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ABSTRACT 

We find that if we live at the center of an inhomogeneity with total density contrast 
\6q\ ~ 0.1 — 0.15, dark energy is not a cosmological constant at 95% confidence level. 
Observational constraints on the equation of state of dark energy, w, depend strongly 
on the local matter density around the observer. We model the local inhomogeneity 
with an exact spherically symmetric solution which features a pressureless matter 
component and a dark-energy fluid with constant equation of state and negligible 
sound speed, that reaches a homogeneous solution at finite radius. We fit this model to 
observations of the local expansion rate, distant supernovae and the cosmic microwave 
background. We conclude that the possible uncertainty from large-scale structure has 
to be taken into account if one wants to progress towards not just precision but also 
accurate cosmology. 

Key words: cosmological parameters - large-scale structure of Universe - cosmology: 
observations 



1 INTRODUCTION 



We liave entered the so-called era of accurate cosmology 
(Peebles 2002). The aim is to understand the composition 
and expansion history of the universe at the percent level. 
In particular, one of the goals is to establish if the observed 
acceleration of the universe (Pcrhnuttcr ct al. 1999; Ricss 
et al. 1998) is driven by the cosmological constant or by 
a fluid with negative pressure, dark energy. At the most 
basic level the question is if the observed equation of state 
w is compatible with —1 or not, the value corresponding to 
the cosmological constant. It is therefore crucial to study 
all possible systematic effects on w (Amendola et al. 2010; 
Sinclair et al. 2010; Marra & Paakkonen 2010; de Lavallaz 
& Fairbairn 2011; Romano & Chen 2011). 

As our observations are confined to the light cone, there 
is an intrinsic degeneracy between temporal evolution and 
spatial variation around us. In particular, inhomogeneities 
around us are degenerate with the properties of dark en- 
ergy, most importantly its equation of state. A clear example 
of how intertwined are attempts to detect any evolution of 
dark energy to large-scale structures is given by the so-called 
"void models" . An observer inside a spherical underdensity 
expanding faster than the background sees indeed appar- 
ent acceleration, thus removing the need for dark energy 
(see e.g. Marra & Notari 2011, and references therein). Void 
models strongly violate the Copernican principle and have 
been ruled out - at least in their simplest incarnation - as 
they predict a too strong kinematic Sunyaev-Zel'Dovich ef- 



fect (Garcia-Bellido & Haugboelle 2008; Zhang & Stebbins 
2011; Zibin & Moss 2011; Zibin 2011). While on one hand 
this strengthens the case for dark energy as the likely ex- 
planation for the acceleration of the universe, on the other 
hand it illustrates how large-scale structure can alter the 
determination of cosmological parameters. Therefore, it is 
necessary to adequately model large-scale structures if one 
has to achieve the grand goal of accurately determining the 
composition of the universe. 

In Valkcnburg (2012b) it was shown by means of mock 
data that a local inhomogeneity, of proportions similar to a 
structure on the surface of last scattering that could cause 
the CMB Cold Spot, can have strong effects on our per- 
ception of the equation of state of dark energy. Here we 
extend that analysis to real data using the model of Marra 
& Paakkonen (2012). More precisely, we consider a wCDM 
model endowed with a local almost-linear inhomogeneity 
surrounding the observer and test it against supernova ob- 
servations, CMB anisotropics and local measurements of the 
Hubble parameter. By wCDM we mean a universe contain- 
ing dark matter and a dark-energy fluid with equation of 
state w. In this way we can show how the observed large- 
scale structure of the universe could impact the reconstruc- 
tion of the dark-energy parameters. 

Following Valkenburg (2012b), we consider an inhomo- 
geneity inspired by the observed Cold Spot in the CMB, 
which has a radius of roughly 5° and a temperature devia- 
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tion of roughly 0(50 ~ 200) fiK} The idea is that the Cold 
Spot is a primary CMB anisotropy due to an object on the 
surface of last scattering, and not a secondary effect caused 
by an object along the line of sight (Tomita 2005; Inouc & 
Silk 2006, 2007; Masina & Notari 2009). Such an inhomo- 
geneity has a radius of roughly 1 Gpc and a density contrast 
today of roughly —0.1 (Valkenburg 2012b). It is therefore too 
shallow to bias the wCDM model to the point of removing 
the need for dark energy, as in the void scenario. Nonethe- 
less, as argued above, such a structure may bias the value of 
the dark-energy parameters to a level that may be important 
if one wants to determine whether dark energy is a cosmo- 
logical constant or not. Moreover, as we briefly argue in the 
body of this paper, structures of radius ~1 Gpc and density 
contrast today ~0.01 ~ which still give an interesting effect 
- are not at more than three times the dispersion of the 
density perturbations arising from a close to scale-invariant 
primordial spectrum. Therefore, the setup considered in this 
paper is not in conflict with standard cosmology and, in par- 
ticular, the Copernican principle (Valkenburg et al. 2012). 

We model the inhomogeneity with a particular case 
of the spherically symmetric solution presented in Marra 
& Paakkonen (2012), which features a pressureless matter 
component and a dark-energy fluid with constant equation 
of state and negligible sound speed Cs- The possibility of 
a dark-energy fluid with negligible sound speed has been 
investigated in the literature under various assumptions. 
This generally requires a non canonical scalar field like k- 
essence and kinetic gravity braiding, as opposed to stan- 
dard quintessence models with canonical scalar fields which 
always have Cs = 1 (see e.g. Creminelli et al. 2009; Bertacca 
et al. 2008; Lim et al. 2010; Bertacca ct al. 2010; Deffayct 
et al. 2010; Li et al. 2011, and references therein). Here we 
choose this particular model because it significantly simpli- 
fies the dynamical equations and the numerical analysis as 
there are no pressure gradients that can generate peculiar 
velocities from an initially comoving motion. We use this 
model phenomenologically so as to minimally extend the 
parameter space of the wCDM model by adding only two 
extra parameters: the radius of the inhomogeneity and its 
overall contrast. All the other initial conditions follow indeed 
rigidly. 

The paper is organized as follows. In Section 2 we go 
briefly through the formalism of the model and its initial 
conditions, and in Section 3 we explain how the cosmological 
data analysis has been performed. We show in Section 4 
that the effect of local inhomogeneity on the dark-energy 
parameters can be important and that the inhomogeneity is 
not unlikely to occur. We conclude in Section 5. 



2 THE MODEL 

We consider the case of an observer located at the cen- 
ter of an inhomogeneous sphere embedded in a fiat wCDM 
universe. The inhomogeneities are given by a pressureless 
matter component and by a dark-energy fluid with constant 



equation of state Wout and negligible sound speed (the sub- 
script "out" refers to values at r > rj, where rt is the co- 
moving radius of the inhomogeneity). To be more precise, 
the sound horizon is much smaller than the inhomogene- 
ity scale considered so that we can set Cs = throughout 
the paper. In terms of the radius dependent non-adiabatic 
equation of state this means that we consider: 



'j{r,t) =- 



pX,out{t) 



Px{r,t) 
that is, the pressure is homogenous: 

px{r,t) = PX,out(t) = Wout PX,out{t) , 



(1) 



(2) 



where r is the coordinate radius, i is cosmic time, p denotes 
energy density, p denotes pressure and the label X refers to 
the dark-energy fluid. Since pressure gradients are absent, 
matter and dark energy evolve along geodesies. Moreover, 
we set initial conditions such that dust and dark energy are 
initially comoving. Therefore, the absence of pressure gra- 
dients implies that peculiar velocities between the two flu- 
ids will never develop and that the matter and dark-energy 
reference frames always coincide. Next, we discuss the equa- 
tions governing the dynamics of the model and the relevant 
initial conditions. 



2.1 Dynamical equations 

We adopt the exact spherically symmetric inhomogeneous 
solution with n perfect fluids presented in Marra & Paakko- 
nen (2012), which we limit to the case discussed above. We 
will now briefly introduce the relevant equations, and we re- 
fer to Marra & Paakkonen (2012) for the general equations 
and more details. 

Using the reference frame of the dust and dark-energy 
components, the metric describing our model is: 



ds^ = -dt^ + X^ll^dr' + Yir,tfdn^ , 



1 — k{r)r^ 



(3) 



where Y{r, t) is the scale function, k{r) is the curvature func- 
tion, AQF' — d9^ + siv? 6 AffP' and we have set c = 1. A prime 
denotes partial derivation with respect to the coordinate ra- 
dius r, whereas a dot denotes partial derivation with respect 
to the coordinate time t. The curvature function is time in- 
dependent because of the adopted reference frame and sound 
speed (Marra & Paakkonen 2012). The metric (3) is writ- 
ten in the same form as the Lemaitre-Tolman-Bondi (LTB) 
metric. However, with the inclusion of a dark-energy fluid, 
the dynamics is no longer that of the LTB metric. The met- 
ric of Eq. (3) reduces to the Friedman-Lemaitre- Robertson- 
Walker (FLRW) metric if k{r) =const and Y{r,t) — ra{t), 
where a{t) is the scale factor. 

We label the dust component with M and the dark- 
energy component with X. The conservation equation for 
the dust source can be solved directly and gives: 



PM{r,t) _ Y''{r,t)Y'{r,i) 
PM(r,t) Y^{r,t)Y'{r,t) 



(4) 



1 See e.g. Cruz ot al. (2006); Zhang 
et al. (2011). 



Huterer (2010); Bennett 



where t is the initial time at which we give the initial con- 
ditions. As explained before, there are no peculiar velocities 
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Figure 1. Shape of the auxiliary function W3 {x,0) that is used 
in Eq. (11) to model the curvature profile. 



between the two fluids and the remaining dynamical equa- 
tions reduce to: 



y'ir,t)^'-m^^kiry 



Y'{r,t)^ 



Y{r,t) 

GF'{r,t) GF{r,t)Y'{r,t) [k{r)r^]' 

Y{r,t)Y{r,t) Y'^{r,t)Y{r,t) 2Y{r,t) 



(5) 



(6) 

F{r,t) ^ - 4-KY''(r,t)Y(r,t) px,out{t), (7) 

px {r, i) = - [px (r, t) + px.out (i)] [Hr (r, t) + 2Ha (r, t)] , 

(8) 

where the radial and angular expansion rates are Hr = 
Y' /Y' and Ha = Y/Y, and F is the total efTective grav- 
itating mass which also satisfies the following consistency 
equation: 



^' (r, t) = 4nY^ (r, t)Y' (r, t) [pM (r, t) + px (r, t) 



(9) 



Eq. (6) is the r-derivative of Eq. (5) and allows us to solve 
directly for the unknown functions Y , Y' , F and px without 
having to take numerical derivatives. If ui = lOout = —1, this 
solution becomes the usual ALTB model which has been 
studied recently in, e.g., Enqvist & Mattsson (2007); Sin- 
clair et al. (2010); Marra & Paakkonen (2010); Valkenburg 
(2012a); Romano & Chen (2011). 

Finally, in this particular case the light-cone equations 
have the same form as in the LTB model: 



dt 
dz 



Y' 



(\^z)Y' 



dr 
dz 



yjl — k{r)r'^ 
{l + z)Y' 



(10) 



2.2 Initial and boundary conditions 



We fix the flat uiCDM background model by setting h, Q,x 
and TOout, where Q.x is the present-day background dark- 
energy density parameter and h is the present-day dimen- 
sionless Hubble rate defined by Hq = 100 ft km s~^ Mpc~^. 
We fix the gauge for the radial coordinate in Eq. (3) 
such that Y(r,t) = r a{t) at some initial time i. We choose 
to parametrize the curvature by means of 



k{r)=KWi[ -,0 



(11) 



where ri, is the comoving radius of the spherical inhomo- 
geneity and: 



Wz {x, 0) 




cos(47rx) 
2 I cos(47r3;) 



+ 



for ^ a: < I 
for I ^ a; < 1, 
for a; > 1 



(12) 

is the third order of the function Wn{x, a), which has been 
defined in Valkenburg (2012a) and interpolates from 1 to 
in the interval a < a: < 1 while remaining C" everywhere. 
Hence k{r) is C^ everywhere, such that the metric is C^ 
and the Riemann curvature is C°. Although the function 
Wz {x, 0) looks rather complicated, it has actually a very 
simple shape as one can see in Fig. 1. The curvature pro- 
file of Eq. (11) is exactly zero for r ^ rt and so the met- 
ric is correctly matched to the exterior spatially-flat toCDM 
model, which means that the central over- or under-density 
is automatically compensated by a surrounding under- or 
over-dense shell. The constant fcc gives the curvature at the 
center of the inhomogeneity and will determine its density 
contrast. A spherical inhomogeneity depends crucially on 
only two physical parameters: the radius and the overall 
density contrast. Therefore, the precise shape of the den- 
sity profile should not be essential and our analysis should 
be representative also of other possible curvature or density 
profiles. 

Next, we have to give initial conditions at t = f for 
F{r,t) = Fm{j') + Fx{r). We choose i such that matter is 
dominant over dark energy, Fx{r) <^ Fj^f{r), and so the 
model becomes the standard dust LTB model. We can then 
link the curvature function k{r) to the initial condition for 
FAiir) by demanding that the universe has the same age i 
for any r. That is, we demand a homogeneous Big Bang, im- 
plying the absence of decaying modes in the matter density 
(Zibin 2008). In particular we can use the following ana- 
lytic result of Van Acoleyen (2008) valid for a linear matter 
density contrast: 



k{r) 



-a^(i)Hout{if 5f,^j{t), 



(13) 



which clearly relates the curvature at the center kc to the 
matter contrast at the center 5m = <5fj,j (0). The latter term 
is the contrast in the gravitating mass and is defined as: 



,ir) 



FM{r) 



Fa 



t(r-) 



1, 



(14) 



where FM,out ~ ^a^{i)r^pM,out{i) is the corresponding 
background gravitating mass and the gauge Y{r,t) = r a{i) 
has been used. The initial matter density is then: 

^m(0 



pM{r,t) 



(15) 



In order to have initial conditions without decaying 
modes in the dark-energy component, we have to set its ini- 
tial profile according to the following relation valid during 
matter domination and c^ ^ 1 (Ballcstcros & Lcsgourgucs 
2010): 



Sfx jr) _ ^ _ 1 + ^"out 

^Fm (^) ^M 1 - SWout ' 

where, analogous to the matter contrast, we have, 

Fx(r) 



Sfx (r) 



Fx,out{r) 



(16) 



(17) 
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Fx,out = ^a^{i)r^px,out{i) and 5x = ^^^(O). The initial 
dark-energy density is then, 

Hence, ah the initial conditions relative to the inhomoge- 
neous patch are indeed specified by a given curvature pro- 
file k{r). 



3 COSMOLOGICAL DATA ANALYSIS 

In this Section we explain how we compare the predictions 
of this model with supernovae, Hubble rate and cosmic mi- 
crowave background observations. We decided not to include 
baryon acoustic oscillations in the analysis as perturbation 
theory in an inhomogeneous background has not been thor- 
oughly understood yet^ (Zibin 2008; Clarkson et al. 2009; 
Alonso et al. 2010). 



3.1 Supernova observations 

We use the Union2 SN Compilation (AmanuUah et al. 2010), 
which consists of 557 type la supernovae in the redshift range 
z = 0.015 — 1.4. As we are considering an almost-linear in- 
homogeneity surrounding the observer, we are not departing 
strongly from the standard model. Therefore it should be a 
good approximation to use the magnitude-redshift and cor- 
relation tables provided by Amaimllah et al. (2010). 



3.2 Cosmic microwave background 

The metric of Eq. (3) is matched to the background FLRW 
metric at a redshift at which radiation is still negligible. In 
this way the last scattering surface, which is responsible for 
most of the CMB anisotropies, is outside the inhomogeneous 
patch and a standard analysis of the primordial CMB power 
spectrum is possible. One has to replace the inhomogeneous 
model with an effective FLRW metric which accounts for 
the different angular diameter distance to the surface of last 
scattering of the CMB as compared to the homogeneous 
background model. This is done by placing an FLRW ob- 
server {So = 0) in the same coordinate system at r = 
but at a different time than to, such that this observer's 
angular diameter distance to the surface of last scattering, 
which lies at some constant time ^ls , agrees with the actual 
LTB observer's angular diameter distance to the surface of 
last scattering. The physics at last scattering itself is unaf- 
fected, since the effective FLRW observer is placed in the 
same FLRW universe in which the LTB patch is embedded, 
albeit at a different time. The CMB spectrum is then cal- 
culated using CAMB (Lewis et al. 2000). See Biswas et al. 
(2010); Moss et al. (2011); Marra & Paakkonen (2010) for 
more explicit details about how the effective model is ob- 
tained. Note that there are other contributions to the CMB 
coming from secondary effects, which in the inhomogeneity 
may differ from those in the effective FLRW metric, and are 
due to the photons traveling through inhomogeneities inside 
the void, such as the late-time ISW effect and weak lensing. 



See, however, Nishikawa et al. (2012) for a recent development. 
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Table 1. Priors imposed on the parameters in the numerical anal- 
ysis. The size of the LTB patch d{ri,) is defined in Eq. (19). The 
additional constraint Qx > with fi^ = 0, in fact implies a non- 
flat prior on both Q^^h'^ and h, as explained in Appendix A. 



For the same reason as for which we ignore the baryon acous- 
tic oscillations, we ignore these secondary effects, since they 
are subdominant and studying them would require knowl- 
edge of the growth of perturbations in an inhomogeneous 
background. We fit the theoretical predictions of our model 
to the WMAP 7-year data release (Komatsu et al. 2011). 



3.3 Parameter estimation 

We perform a Markov-Chain Monte-Carlo likelihood anal- 
ysis using COSMOMC (Lewis & Bridle 2002). We calculate 
all distance measures using an improved version of VoiD- 
DlSTANCESlI^ (Biswas et al. 2010), which now acts as a 
wrapper around CAMB (Lewis ct al. 2000), necessitating no 
changes to CAMb's source code and minimal changes to Cos- 
MOMC's source code. We combine this module with the 
ALTB module ColLambda* (Valkenburg 2012a) for calcu- 
lating all metric functions, which we extended to include 
the numerical solutions to the scenario discussed here, with 
Wout 7^—1 and Cs — 0. With this setup, for every selected 
vector of parameter values, we calculate the theoretical pre- 
dictions for supernova distances, the local Hubble rate and 
the CMB power spectrum. For reference, next to the inho- 
mogeneous model we analyze its homogeneous background 
model, wCDM, which is described by the same model but 
has 26 = and So = 0. 

We take fiat priors on the parameters listed in Ta- 
ble 1. Most of these are the usual cosmological parameters: 
the background present-day Hubble rate h, the background 
baryon density flhh^, the background dark-matter density 
^dmh^ , the equation of state of dark energy lUout, the op- 
tical depth to re-ionization r, the amplitude of primordial 



^ http : //web . phys ik . rwth- aachen . de/download/valkenburg/ 
■* http : //web . physlk . rwth- aachen . de/download/valkenburg/ 
ColLambda/ 
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scalar perturbations As, the tilt of the spectrum of primor- 
dial scalar perturbations ns and its running as- We set the 
spatial curvature outside the LTB patch, Qt, as well as the 
amplitude of primordial tensor perturbations to zero. 

As discussed in Section 2.2, two additional parameters 
describe the LTB patch: the curvature at the center kc and 
the comoving radius r^. As shown in Table 1, we will use 
as actual parameters the present-day total density contrast 
and the present-day proper size of the radius, respectively. 
The latter is given by 

'■'■^ Y'{r,to) 



din) 



y/I — k{r)r^ 



dr , 



and we define the former as 

Pi. 



5o = 



Pout 



(19) 



(20) 



max(pin,Pout) 

where pin = Pm(0, to) + Px{0,to) and pout = PM,out(to) + 
pj!f,out(io)- We calculate the contrast using the total den- 
sity as in principle the dark-energy component can be as 
inhomogeneous as the matter component, and So should be 
a relevant physical quantity to be used. We chose a some- 
what unusual definition by using max(pin, Pout) in the de- 
nominator. With this definition, this quantity fundamen- 
tally satisfies —1 < So < 1, such that the parameter- 
space volume in over- and under-densities is equally dis- 
tributed. With the more usual definition we would have 
had — 1 < P'°~P°"* < (X), which in the Bayesian parame- 
ter estimation induces a strong prior favouring large over- 
densities, possibly excluding under-densities from the anal- 
ysis. In practice, however, the preferred values for this pa- 
rameter are small, such that the difference between both 
definitions is almost negligible. 

We discuss the priors in more detail in Appendix A. 
In the next section, where we discuss the results of the 
MCMC parameter estimation, we explore the full param- 
eter space, never fixing the background parameters to some 
central value. Therefore we can always marginalize over all 
parameters, and do not bias the result in any way. That is, 
any possible degeneracy, expected or unexpected, between 
the cosmological parameters and the LTB parameters will 
show up and will not influence the results without being 
noticed. 



4 RESULTS 

The four left panes in Figure 2 show the parameters on which 
the main focus in this paper lies. This figure shows the two- 
dimensional marginalized posterior probability distributions 
of the background dark-energy density Six, the dark-energy 
equation of state Wout, the boundary of the inhomogene- 
ity in redshift space zt, and the total density contrast in 
the inhomogeneity So (see Eq. (20)). Alternatively, in the 
four right panes of Figure 2 we use as proxies for zt and 
So the apparent size that the inhomogeneity would subtend 
if located at the last scattering surface and its correspond- 
ing temperature anisotropy, respectively. Here we define the 
temperature perturbation AT/T as the relative difference in 
CMB temperature at the center of the inhomogeneity and 
outside, in the homogeneous (average) background. This is 
not necessarily representative for the average temperature 
in the spot. 



The most interesting result is the clear degeneracy be- 
tween So and Wout in the lower right pane in the left of Fig- 
ure 2, and the degeneracy between So and fix in the upper 
right pane. This graph explicates the necessity of properly 
modeling the inhomogeneity of the local universe, before any 
conclusion can be drawn on the properties of dark energy, 
in particular about the fundamental value of its equation of 
state: if the local density is ignored, the value of Wout can be 
misestimated by possibly 50%. Note also that the redshift 
up to which the inhomogeneity extends, zt, is hardly of any 
infiuence on the central value of Wout or Qx ■ The bias on the 
parameters is indeed of opposite sign for opposite So- There- 
fore, as we marginalize over So in the combined posterior of 
Wout and Zb (lower left pane in Fig. 2), the total effect is com- 
pensated (we will come back to this point with Figure 3). 
Note, however, that the scatter does increase with the size of 
the inhomogeneity. The right of Figure 2 shows the dimen- 
sions that the inhomogeneity would have on the observed 
CMB temperature map, if it were centered on the observer's 
surface of last scattering. In this situation there are hence 
two identical inhomogeneities in the universe: one surround- 
ing the observer, one centered on the observer's surface of 
last scattering. This figure shows that even spots that do 
not violate the observables of Section 3, do induce a strong 
bias on Wout, following the findings of Valkenburg (2012b). 
As before we see that the apparent size of the inhomogeneity 
has almost no effect on the value of Wout or Qx - 

In Figure 3 we show the one-dimensional marginalized 
posterior probabilities for fix and uiout, under different pri- 
ors on So, imposed by means of importance sampling on 
the MCMC chains that explore the full range of So- The 
result is displayed in black, and for comparison the con- 
straints on these parameters for the homogeneous luCDM 
model {zb = 0, So = 0) are displayed in dashed red lines. The 
top row shows Jlx and Wout marginalized over all values of 
So, both positive and negative. Since positive and negative 
values of So have opposing effects on Wout and Qx , marginal- 
izing over So mostly widens the tails of the distributions. 
This is already an important result showing how inhomo- 
geneities contribute to the error budget in the cosmological 
parameters. This uncertainty from large-scale structures is 
expected to become more important when future data will 
tighten the confidence regions of the parameters of interest. 
If we impose, however, the prior that So > (second row 
in Figure 3) or Jo < (third row from top), we find even 
stronger results: we see indeed that both the tails and the 
central values of Q-x and Wout shift. This shift is significant 
if one wants to progress towards not just precision but also 
accurate cosmology. If we push the magnitude even further, 
pretending we know that the local density must be either 
So =0.1 or So = —0.15, as in the bottom row in Figure 3, 
then we find that w; = — 1 is excluded at 95% confidence 
level (c.l.) in both cases: w < —1.03 and w > —0.98, respec- 
tively. However, these particular models are only included 
at 99.7% c.l. 

In order to better show the degeneracy of Qx and Wout 
with So we plot again in Figure 4 the corresponding two- 
dimensional marginalized posterior probabilities with 68%, 
95%, 99.7 % and 99.99% confidence level contours. Also plot- 
ted for comparison are the 95% c.l. one-dimensional con- 
straints on Qx and Wout for the standard wCDM model. 
This plot is meant to justify the values So — 0.1 and 
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Figure 2. Two-dimensional marginalized posterior probability distributions for the most interesting parameters characterizing the tuCDM 
model endowed with a local inhomogeneity considered in this paper. The parameters are constrained by the Union2 SN Compilation 
(AmanuUah et al. 2010), the WMAP 7-year CMB power spectrum (Komatsu et al. 2011) and the recent determination of the Hubble rate 
by Riess et al. (2009). Inner tangerine-tango coloured contours are 68% confidence level (c.l.) contours, while the outer dark-red coloured 
contours are 95% c.l. The left and right figures show the same information, however the redshift of the boundary of the inhomogeneity, 
2;,, and the total density contrast, (5o, on the left, are on the right traded in for the angular diameter that the inhomogeneity would have 
if it were located at the observer's last scattering surface and the temperature fluctuation that it would induce, respectively. 



5o = —0.15 used in the bottom row of Figure 3. It shows 
indeed that large values of |(5o| that can significantly bias 
fix and lOout are still within the 99.7% c.l. 

4.1 Sensitivity to local Hubble-rate constraints 

In Figure 5 we compare the effect of different Hubble-rate 
observations on the resulting posterior probabilities of Q.x 
and Wout, when we fit the luCDM model endowed with a lo- 
cal inhomogeneity to Hioc, SNe and CMB. As explained in 
Section ??, we compare the three different values from Riess 
et al. (2009) (solid black), Freedman ct al. (2001) (dashed 
red) and Sandage et al. (2006) (dashed dotted blue). We 
see that the constraints on fix do depend on the chosen 
measurement for -ffioc, while the resulting constraints on 
uiout hardly depend on -ffioc. This should be expected as 
SNe observations (constraining both fix and Wout) are in- 
sensitive to H\ac while CMB observations (constraining Q.x 
but weakly uiout) are instead sensitive to Hioc. These results 
imply that our conclusions with regard to Wout are robust 
against different observational determinations of i/ioc, and 
the constraints hence mostly follow from the SNe and CMB. 

4.2 kSZ effect 

In Figure 6 we show the CMB dipole that observers at differ- 
ent radii would observe, for a given configuration; an over- 
density on the left, an under-density on the right, in both 
cases the LTB patch has a radius of 1 Gpc, but with differ- 
ent values for lOout. Both cases fit the data roughly as well 
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Figure 5. One-dimensional marginalized posterior probabilities 
of Qx and «)out for the inhomogeneous model, given CMB, SN 
and iJioc observations, comparing different constraints on -ffioc- 
Riess et al. (2009) (solid black), Freedman et al. (2001) (dashed 
red) and Sandage et al. (2006) (dashed dotted blue). The conclu- 
sions about the effect of the inhomogeneity on uiout are robust 
against different observational determinations of -ffioc- 



as the standard wCDM, while still giving an interesting bias 
on lUout. 

We obtained these figures by - at each radius - start- 
ing an integration of the geodesic equations in two direc- 
tions (negative and positive r-direction), back to the surface 
of last scattering, which lies at constant time in the syn- 
chronous gauge of the LTB metric. The difference in redshift 
to the surface of last scattering in both directions is then 
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Figure 3. One-dimensional marginalized posterior probabilities for Qx and Wout, constrained by CMB, SNe, and -ffioc- The standard 
uiCDM constraints on these parameters are shown in red dashed lines and are the same in the different panes. The blue vertical line 
serves as a guide for the eye, always going through the maximum of the uiCDM value. The constraints on the uiCDM model endowed 
with a local inhomogeneity considered in this paper are given under different priors on Sq: —0.2 < So < 0.2 (top), < <5o < 0.2 (second 
row from top), —0.2 < <5o < (third row from top), and for two constraining priors on Sq (bottom): Sq = 0.1 (solid black line) and 
So = —0.15 (dashed-dotted black line). A prior on So that averages around zero, widens but does not shift the posterior distributions. If 
we know instead the sign of 5o, constraints on tUout and Qx shift by as much as 5% ~ 10%. In the bottom row, for So = 0.1 the 95% 
c.l. upper bound on tOout is -1.03. For So = 0.15 the 95% c.l. lower bound on uiout is -0.98. Both priors hence rule out the cosmological 
constant at 95% c.l., given current observations. 



translated into a ATcMB/rcMB = {z+ — Z-)/{2 + z+ + Z-). 
Only for small radii this is to a good approximation equal to 
the dipole in a spherical harmonics expansion of the CMB 
temperature map. On the vertical axis on the righthand 
side we list the corresponding peculiar velocity that an ob- 
served temperature difference corresponds to, if it were the 
effect of peculiar velocity alone. For both panes, left and 
right, the peculiar velocities do not exceed the magnitude 
of expected random peculiar velocities. Therefore the kine- 
matic Sunyaev-Zel'dovich effect that such velocities induce 
on CMB photons (Garcia-Belhdo & Haugboelle 2008; Zhang 
& Stebbins 2011; Zibin & Moss 2011) should be at present 
undetectable. So as to strengthen this claim it is useful to 
look at the findings of Valkenburg et al. (2012) where con- 
straints on the ALTB model from kSZ observations have 
been computed. While in the present paper dark energy is 
not the cosmological constant, the above analysis should give 
nevertheless an estimate of the kSZ signal. Therefore, this 
suggests that structures with a contrast of roughly ~0.1 ex- 
tending for a radius of 1-2 Gpc are not excluded by present 



observations. A thorough study of the kSZ effect in these 
models is left to future work. 



4.3 FLRW Observer's wiz) 

Following Clarkson et al. (2007) one can, given a luminosity 
distance-redshift relation in a homogeneous universe, com- 
pute what the underlying w{z) of the dark-energy fluid is. In 
the homogeneous universe (described by the FLRW metric) , 
one can find indeed an exact relation between w{z) and the 
first and second derivatives of the luminosity distance with 
respect to redshift and two more parameters, Qk and fim. 
If an observer knows the latter two parameters from other 
observations, and deduces the first and second derivatives of 
the luminosity distance from SN observations, the observer 
can derive ^1(2). In the scenario studied here, at background 
level 10 is not a function of time or redshift. However, the 
inhomogeneity comes into play in the luminosity distance- 
redshift relation. Therefore, an observer that falsely assumes 
that the metric surrounding him/her is FLRW will in fact 
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towards values of «i that are far from — 1 . 



9.1-10 



-9.1-10 



5o = 0.017, Wo,t = -1.08 




-900 -450 450 
d(r) [Mpc] 



900 



5o = -0.030, WoLit = -0.97 



1.7-10 




-1.7-10 



-900 -450 450 
d(r) [Mpc] 



900 



Figure 6. The CMB-dipolc observed by observers at different radii d{r) in an LTB patch with a radius of 1 Gpc, for an over-density 
(left) and an under-density (right). On the right vertical axis we list the corresponding peculiar velocity that is derived assuming that 
the observed CMB dipole is caused solely by the peculiar velocity of the observer. The magnitude of the velocities does not exceed the 
magnitude of expected random velocities. See Section 4.2 for more details. 



see a redshift dependence in w. We calculate the observed 
Wobs{z) (Eq. (3) in Clarkson ct al. 2007) for the two exam- 
ple models of Fig. 6, and show the result in Fig. 7. Inside 
the inhomogeneity, Wohsiz) shows a very clear signature of 
the matter distribution: the contracting core and expanding 
compensating shell for the over-density show corresponding 
effects on uiobsC^). The inverse holds for the under-density. 

Therefore, if one performs an analysis such as in 
Shafieloo et al. (2009); Zhao ct al. (2012), one may find 
a significant deviation from a constant w, while fundamen- 
tally w is constant at the background level. In particular, it 
is very interesting to note that the w{z) reconstruction by 
Zhao et al. 2012 (see e.g. the pane (A2) of Fig. 1 in that 
reference), if interpreted within this framework, could indi- 
cate the presence of a large-scale underdensity around us. 



and not of a possibly time-dependent equation of state. It 
is indeed worth noting the similarity of the result by Zhao 
et al. (2012) with the observed Wobs(z) shown in Fig. 8, 
which corresponds to a ACDM model endowed with a lo- 
cal underdensity of central contrast So = —0.06 and redshift 
boundary Zb = 0.4. 



4.4 Other parameters 

In Figure 9 we show posterior probabilities of Ho, which is 
the expansion rate of the background universe (and does 
not correspond to the locally observed expansion rate), 



So 



max(pi„,pout) ' 



1, and Sx,o 



All other parameters that we allowed to vary, listed in Ta- 
ble 1, show practically no deviation in their constraints in 
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Figure 7. The function w^hai^) eis defined in Clarkson ct al. (2007), wliich is the equation of state of dark energy that an observer thinks 
to see if the observer falsely assumes that the universe is described by the FLRW metric. The examples shown here are the same two 
cosmologies as in Fig 6; an over-density (left) and an under-density (right). In both cases the observed w^hai^) (solid red) matches closely 
with the fundamental uiout (dashed blue) once the radius is reached where the spherically symmetric metric matches to the surrounding 
FLRW metric. The inhomogeneity causes clear features in «)obs(^). the contracting core of the over-density increases the magnitude of 
Wohai^)t while the under-dense compensating shell has the opposite effect. For the under-dense center, the inverse holds. In this picture 
u>obs(^) is computed using the exact solutions, while taking second derivatives of observed distances would not necessarily reveal these 
features. 



Comparison to Zhao et al. (2012) 
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Figure 8. As in Fig. 7 but for a ACDM model endowed with 
a local underdensity of central contrast 5q = —0.06 and redshift 
boundary Zf, = 0.4 such that the observed i«oba(^) appears qual- 
itatively similar to the w{z) reconstruction by Zhao ct al. (2012) 
(see e.g. the pane (A2) of Fig. 1 in that reference). If interpreted 
within this framework, the results of Zhao ct al. (2012) could indi- 
cate the presence of a large-scale underdensity around us, rather 
than a possibly time-dependent dark-energy equation of state. 



the presence and absence of the inhomogeneity. The con- 
straint on Ho is significantly weakened when one takes into 
account the possibility that we may live in a local inhomo- 
geneity. This result is similar to the findings in Valkcnburg 
& Bjaelde (2012), even if the local inhomogeneity consid- 
ered is orders of magnitude different. Moreover, the actual 
matter perturbation that is allowed by the data can be as 
large as \6m,o\ — 0.5, leaving the total density perturbation 
around \So\ ~ 0.1, because the energy perturbation in the 
dark-energy fluid is generally small as Wout is never very far 
from -1. 



4.5 Probability under homogeneous initial 
conditions 

Given the fact that we observe several large spots in the 
CMB, and the possibility that these spots are the result 
of density perturbations on the surface of last scattering, 
we can argue that living in such a perturbation must have 
a non-zero probability.'' If all perturbations arise from a 
smooth, close to scale invariant spectrum of primordial per- 
turbations, it is not obvious that large cold and hot spots 
should exist, and there is an ongoing debate about this topic 
(Cruz et al. 2006; Zhang & Hutcrcr 2010; Ayaita ct al. 2010). 
Let us nonetheless quantify the probability of having 
the density perturbations that we took as examples for 
Figs. 6 and 7. If these perturbations come from the same 
spectrum as the perturbations that we observe in the CMB, 
then the probability of their existence can be approximated 
by the variance of the gaussian density field, smoothed by 
a top hat filter with a radius that corresponds to the ra- 
dius of the density perturbation under consideration (Kolb 
& Turner 1990). It must be noted that, because of the com- 
pensated shape of the density profile that we consider, tak- 
ing the full radius of the spherical patch, rt, as the radius 
of the top hat filter would give almost zero density pertur- 
bation. Therefore we choose the radius at which the den- 
sity changes sign as the smoothing radius. This is roughly 
at rb/2, but we use the numerically obtained exact value. 
Comparing at the time of decoupling when dark energy is 



^ We would like to point out that for the almost-linear models 
considered in this paper the observer does not need be very close 
to the center so as not to see a too large CMB dipole. For the cases 
of Fig. 6, for example, the dipole is never larger than the observed 
value of ^ 10"''. For larger contrasts there will be regions where 
observers would sec a larger-than-measured dipole, but these re- 
gions will not occupy the majority of the inhomogeneity. Finally, 
we would like to stress that we have placed the observer at the 
center simply to simplify the numerical calculations. 



10 Marra, Pddkkonen and Valkenburg 





-0.2-0.1 0.1 



-0.5 





-0.04-0.02 



Figure 9. Onc-dimcnsional marginalized posterior probabilities of Hq, which describes the age of the universe and does not correspond 
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structure weakens significantly the bounds on the background expansion rate (in red we show the constraint on Hq in the homogeneous 
uiCDM model) . Secondly, the local energy perturbation consists of an almost negligible dark-energy perturbation and a significant dust 
perturbation, which is not apparent if one considers the total Sq alone. 



negligible, and with the primordial spectrum of perturba- 
tions of the two models respectively (since the models were 
fit to the CMB, they carry spectral parameters), we find 
that the over-dense model considered in Figs. 6 and 7 is at 
three times the dispersion of the smoothed density field of 
its cosmology (CMB spectrum), and the under-dense model 
is at six times the dispersion of its cosmology. Notably the 
over-dense model is not at all unlikely to occur, while it does 
give a large effect on both Wout and Wohs as shown in Fig. 7. 



5 CONCLUSION 

We have analyzed present observations of the local expan- 
sion rate, distant supernovae and the cosmic microwave 
background within a flat wCDM model endowed with a local 
almost-linear inhomogeneity surrounding the observer. We 
have found a siginficant impact on the dark-energy parame- 
ters, in particular on the equation of state which is strongly 
degenerate with the inhomogeneity contrast. The implica- 
tions of this degeneracy are twofold. On one hand we have 
shown that with prior knowledge on the inhomogeneity, to 
be obtained possibly with some future probe, it is already 
possible to rule out the case of the cosmological constant 
with current data. On the other hand, even if future probes 
exclude the case of the cosmological constant in a homo- 
geneous universe, this still may be due to a poor modeling 
of the large-scale structure of the universe. The same con- 
clusions apply to constraints on the time variation of the 
equation of state. 

The analysis in the present paper is but a first step to- 
wards a more accurate reconstruction of the cosmological 
parameters. We have indeed chosen, for technical reasons, 
a very specific dark-energy model and inhomogeneity pro- 
file. Before drawing definitive conclusions, a more compre- 
hensive analysis should be performed. Firstly, it would be 
particularly interesting (even though perhaps challenging) 
to consider a nonzero dark-energy sound speed. However, 
since we found that the dark-energy component is only very 
mildly inhomogeneous, we do not expect our results to be 



strongly dependent on the assumption of a negligible sound 
speed. Secondly, it would be interesting to consider more 
inhomogeneous patches with more general density profiles, 
possibly with the observer at randomized positions (Marra 
ct al. 2007, 2008; Valkenburg 2009; Szybka 2011; Flanagan 
et al. 2012). 
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APPENDIX A: EFFECT OF CONSTRAINTS 
ON PRIORS 

In Figure Al we show the prior probability of ftdmh^, Ho 
and derived parameters, since they dependent non-linearly 
on combinations of input parameters on which we take flat 
priors, and hence their prior probability is non-flat. We ob- 
tained these prior probabilities by running the MCMC anal- 
ysis without any data, accepting all points in parameter 
space with equal probability. Maybe surprisingly, the ac- 
tual prior probabilities on r^dm'*^ and Hq are not flat, even 
though we do list their prior ranges in Table 1 as flat, and 
we indeed gave flat priors on these parameters in the in- 
put of the MCMC simulation. The non-flatness stems from 
the additional constraint fix > 0, demanding that the en- 
ergy density of the dark-energy fluid is positive (if the Dark 
Energy is a pure cosmoillogical constant nothing prevents 
f^A < from happening). Because we set flk ~ 0, this con- 
dition is satisfied only when flm < 1, or flmh^ < h? ■ 

To explain the relation, let us simplify the priors to 
< Q,rah? < 1 and Q < h < 1. The starting point is 
a flat prior on these parameters, i.e. that the probability 
P(f2ni'i^) = constant and P{h) = constant, normalized such 
that J P{x)dx = 1 for x being both h and Q,yah? . If we 
write shorthand notation C for the condition that Q,x > 
and Q,k ~ 0, then imposing C we have the conditional joint 
probability. 



p(a„ft^;^)oc| ^'("-'^''^l^)^ 


= constant 


if C is true 
if C is false ' 

(Al) 


normalized such that 






rl /■! 







P{nn,h\ h\C) d{n,Ji- ) dh = l. (A2) 



Then we find for the probabilities of h and fidm'i^, up to 
normalization constants, 

P{h\C)= I P{n^h^,h\C)d{nmh^) 
Jo 



ex / d(f2mft ) (X h , 
Jo 

P{Qmh^\C) = [ P{n„,h^,h\C)dhoc f 

Jo Jh=^n,„h'- 



ex 1 — V Slnift^ ; 



(A3) 
dh 

(A4) 
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Figure Al. Priors on fi^m'^^i ^0 a-nd derived parameters, displaying tliose tiiat are non-flat. The black (solid) lines represent the 
priors when the inhomogeneous metric is used, while the blue (dash-dot) lines represent the priors in the homogeneous uiCDM case. 
For comparison, in red (dashed) we show the constraints on Q^^h^ and Hq when the homogeneous toCDM case is fit to the data. The 
narrowness of the red curves shows that the data are very constraining on these parameters, and so the non-flat prior is of no importance. 
These prior distributions are obtained by running the MCMC analysis without any data, that is, accepting all points in parameter space 
with equal probability. 



in agreement with the priors in Figure Al, favouring large 
h and small Q^nh^ ■ However, as can be seen from the red 
dashed curves in Figure Al, the data are constraining Qmh^ 
and Ho so tightly, that the non-flatness of the prior has no 
effect on the final parameter estimation. 
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